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Abstract. 

A recent advance in improving the spatial resolution of magnetic force microscopy 
q \ (MFM) uses as sensor tips carbon nanotubes grown at the apex of conventional silicon 

cantilever pyramids and coated with a thin ferromagnetic layer [1] . Magnetic images of 
high density vertically recorded media using these tips exhibit a doubling of the spatial 
frequency under some conditions [T] . Here we demonstrate that this spatial frequency 
doubling is due to the switching of the moment direction of the nanotube tip. This 
\ results in a signal which is proportional to the absolute value of the signal normally 

' observed in MFM. Our modeling indicates that a significant fraction of the tip volume 

t^J- . is involved in the observed switching, and that it should be possible to image very high 

' bit densities with nanotube magnetic force sensors. 

i> " 
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* * 

PACS numbers: 75.75.+a,78.67.Ch 
Spatial period doubling has been observed for several carbon nanotube tips with 
different track widths and bit densities. The MFM images reported here were made 
using the nanotube tip shown in the inset of Figure [U It was approximately 250nm long, 
with a ferromagnetic coating to a total tip diameter of 16nm. The means of producing 
metal-coated carbon nanotube tips on AFM cantilevers and the techniques for imaging 
magnetic media using these tips have been described previously [U [2] . Briefly, carbon 
nanotubes were grown using wafer-scale chemical vapor deposition at the apexes of the 
pyramids of commercial silicon tips intended for tapping mode atomic force microscopy. 
For the present measurements the nanotubes were shortened to a length of about 250 
nm using an electrical cutting method [3], aligned approximately perpendicular to the 
cantilever using a focused ion beam [2], and then coated to a total thickness of about 
16 nm with a Ti/Co/Ti trilayer by e-beam evaporation from a direction parallel to the 
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nanotube long axis (see the inset in Figured]). The magnetic imaging was done using the 
Tapping/ Lift™ mode of a Digital Instruments Nanoscope III SPM at room temperature 
in air. In this mode, the topography of the sample is first determined for each line with 
a scan at low tip-sample spacing z , then the tip is retracted a specified distance and 
a second line scan is made while recording the deviation of the phase angle S of the 
cantilever response with the cantilever driven slightly below its resonance frequency. 5 
is proportional to dF z /dz , the derivative of the force on the cantilever with respect to 
z . 

The images presented here were made on a vertically polarized magnetic medium. 
A collection of phase shift images at 300 kilo-flux changes per inch (kfci) and different 
tip-sample spacings (zo's) are shown in Figure [2(a), and at constant lift height (15nm) 
and different bit densities in Figure [2](b) . Cross-sections of the data through the center 
of the tracks are displayed in Figure E](c,d). Modeling of this data as described below 
is shown in Figure E](e,f). At high values of z and high bit densities the phase images 
have the same period as written, but at low z 's and low bit densities anomalies in 
the images gradually develop into sharp features with double the original periodicity. 
These anomalies are due to switching of the orientation of the magnetic moment of 
the tip. This can be demonstrated most convincingly by inspecting the images and 
cross-sections of Figure [2](b,d). In this case a background from the sections of the image 
without written bits has been subtracted out, and it can be seen that at low bit densities 
the phase shift always stays below the average background level, keeping the tip-sample 
force attractive by switching of the tip moment direction just as the force derivative 
crosses zero. 

This conclusion is supported by detailed modeling: We assume that the magnetic 
medium is composed of slabs of length s in the x direction, height h in the z 
direction, and width W in the y direction (Figure DD(a)). The slabs have a uniform 
magnetization with moment direction alternating between parallel and anti-parallel to 
the z-axis direction. The tip is assumed to have a square cross-section with width w 
and length L, with the tip end a distance z from the upper surface of the medium. 
The magnetic fields above the sample (displayed as field lines in Figure 11(b)) were 
calculated both analytically and numerically. In our numerical modeling the vector 
between the individual medium dipole moments [x m ,y m , z m ] and a position [x, y,z] is 
f = (x — x m )x + (y — y m )y + (z — z m )z. Then the z and x-component of the field from 
the individual dipoles is given by 

D M m m / '3(z - z m f \ 

B *& = 1^3" [— 1J 

with r = \r\. B y is given by setting x —> y and x m —> y m in the second equation. 
If we assume that the tip can also be represented by a sum of point dipoles m t = m t h, 
where ft is a unit vector in the tip moment direction, the force gradient on the tip is 
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Figure 1. (a) Tip and sample geometry in model. The inset is a scanning electron 
microscope image of the cobalt-coated nanotube tip. (b) Calculated field lines for a 
300 kfci track above the center of the domains in the magnetic medium, using the 
parameters described in the text. The arrows indicate calculated moment orientations 
for a series of tip-sample spacings of 15, 22.5, 28, 33, 39, and 45 nm using the model and 
parameters described in the text. At small tip-sample spacings the tip moment flips 
to be always anti-aligned with the moment directly below the tip. At large tip-sample 
spacings the tip moments are only slightly perturbed by the medium fields. 
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where r tm is the distance between the individual dipole elements in the tip and medium, 
and (f) is the angle between the tip moment direction and the z-axis. 
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Figure 2. (a) Collection of phase shift images for a 300 kfci track of alternating 
vertically recorded media moments at different lift heights. Each line of the image has 
an offset such that the average phase shift value for that line is zero, (b) Image of a 
section of vertically recorded media with a series of tracks with different bit densities, 
at a lift height of 15 nm. A background has been subtracted from this image such 
that the regions between the tracks have zero averaged phase shifts, (c) Cross-sections 
through the centers of the tracks in (a), (d) Cross-sections through the centers of the 
tracks in (b). The zero phase shift levels for each cross-section are indicated by dashed 
lines, (e) Modeling of (c) as described in the text. The lines in (c) and (e) have been 
offset vertically for clarity, (f) Modeling of the tracks in (d) as described in the text, 
with the zero phase shift levels indicated by dashed lines. 
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For positions near the center of the tracks we obtained analytical expressions for 
the magnetic fields by assuming that the magnetically oriented slabs have infinite width 
W in the y direction. If we take the boundary condition 

H z (r)\ z=0+ - H z (r)\ z=0 - = H z (r)\ z= _ h - - H z (r)\ z= _ h+ = a(r), (3) 

where a(f) is the surface magnetic charge, the magnetic field above the sample can be 
written as [4] 



H z (r,z) 



4tt 2 



where 



Af z = e~ kz (l - e~ kh ) / d 2 ra{f)e ik - r . 

J — OO 



(4) 



(5) 



Taking the surface magnetic charge a(f) to be uniform in y, 

a(r) = M 2ns < x < (2n + l)s 
= - M (2n - l)s < x < 2ns, 



(6) 



We find 
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dQ/dz, $ a scalar potential, the x-component of the field can be written as 
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(7) 



Since H 7 



H x = dQ/dx, which leads to 
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H y is zero by symmetry. Figure [3] displays the calculated fields in the z and x directions 
(a,c), the tip moment orientation angle (j) (b), and the switching fields < H z > c and 
< H x > c (c) for the best fit to the data as described below. It is interesting to note 
that the discontinuous switches of (ft are controlled predominantly by the size of the x 
component of the field. This can be understood by examining the trajectory in field 
that the tip takes in moving from one domain to the next. The ovals in Figure [3b are the 
calculated trajectories of < H z > vs. < H x > for the tip heights listed in the caption. 
The "asteroid" is the critical field calculated for the energy functional form of Eq. 1 
of the main text and has the form first predicted by Stoner and Wohlfarth[5j. The tip 
moment direction is predicted to switch when the field trajectories cross the critical field 
asteroid (solid dots in Figure [3b) . This happens for relatively large values of | < H x > \ 
and small values of | < H z > |. When the field trajectory crosses the Stoner- Wohlfarth 
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Figure 3. (color online) (a) Calculated magnetic fields, averaged over the tip volume 
and multiplied by the tip saturation magnetization M t divided by the anisotropy 
parameter K±, above a vertically recorded magnetic medium with infinite extent in 
the ^/-direction, 12 nm thick in the z-direction, with magnetization oriented in the 
^-direction and alternating in the x-direction with period 2s = 185 nm, for spacings 
between the bottom of the tip and the medium of zo = 15, 22.5, 26, 28, 30, 33, 39, and 
45 nm. The tip has an assumed square cross-section 16nm on a side, and the averaging 
is over a tip length L\ = 48nm. The calculated magnetic fields < H z > normal to 
the magnetic medium are offset by 3 units. The best fit values M t M m / K\ = 17.5, 
<j w / hqL\Ki = (see main text Figure 3) were used for these calculations, (b) 
Calculated variation of the tip moment orientation angle <\> relative to the z axis, 
(c) The ovals show the calculated trajectories of < H z > M t /K\ vs < H x > M t /K\ 
for the various tip heights. The "asteroid" plots the values of the critical fields. The 
intersections of these two sets of curves, indicated by solid symbols, are the fields at 
which switching occurs in the simulations. 
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asteroid at large values of | < H z > | the tip has already switched to the low energy 
configuration. 

For numerical work the medium was assumed to be composed of a collection of 
individual dipole moments fa m m z, where m m ±M m v, with M m the saturation 
magnetization and v the volume of the individual medium elements. In what follows 
we take both the tip and medium volume elements to be cubes 4 nm on a side. This 
results in agreement to within a few percent between our numerical work and analytical 
expressions for the medium magnetic fields. Halving the size of the volume elements (to 
cubes 2 nm on a side) changes the calculated force derivative curve for z = 30 nm (see 
Figure Eb) by about 2%. 

To model the dynamics of the tip flip process, we conceptually divide the tip into 
two domains, one with length L\ close to the medium, the other with length L — L\ 
further away. Each has sufficiently strong exchange fields that the entire volume within 
each domain has the same moment orient ation[6j. The section of the tip furthest from 
the medium is assumed to have its moment parallel to the z axis; that closest to the 
medium has its moment at an angle cj) relative to the z-axis (Figured]). Then the energy 
of the tip in an external magnetic field can be written as: 

-(1 — cos 0) + sin 2 (ft 



M t 



HoLxKx 2 

(9) 



+ — (< H z > cos (fi+ < H x > sin 4 
Ki 

where we have taken the simplest non-trivial forms for the domain wall energy (first 
term) and the anisotropy energy (second term)[5l [7j. The third term in Eq. [9] is the 
energy of the dipole moments of the tip in the external magnetic field. Here [IqKi is 
the anisotropy energy density, a w is the domain wall energy per unit area, M t is the tip 
saturation magnetization, and < H z > and < H x > are the magnetic fields in the z and 
x directions respectively averaged over the tip volume from z — zq to z — zq + L\. To 
simulate the magnetic force images, the tip moment is at first assumed to be parallel 
to the z-axis. The tip is moved to a new position, the local fields are calculated and 
averaged over the tip volume, (j) is moved to the new local minimum in energy (Eq. 
[9j), the force gradient is calculated, and the process is repeated. This modeling results 
in the cross-sections displayed in Figure E](c,f), which reproduce the absence of tip 
switching at high bit densities and high lift heights, and the presence of tip switching 
at low bit densities and low lift heights. When tip switching occurs, the modeling also 
reproduces the fact that the tip-sample force gradient always stays negative, with the 
tip moment reversing as the z-component of the field crosses zero. The quantitative 
interpretation of MFM images in the presence of tip switching is straightforward once it 
is recognized that the phase shift 8 is proportional to the negative of the absolute value 
of the tip sample force gradient (-\dF z /dz \). The experimental phase shift oscillation 
amplitudes decrease much more rapidly than the modeling for bit densities above about 
500 kfci (Figure EF). We believe that this is because the as written bits do not have 
as abrupt moment orientation reversals as our idealized model. Our modeling indicates 



Moment switching in nanotube magnetic force probes 



8 



that nanotube tips with the geometry of Figured] could be used to image bits with sharp 
moment direction transitions with densities above 2000 kfci (13 nm/flux reversal). 

Figure [4] compares the maximum minus the minimum value for S along a cross- 
section through the center of the bits at a bit density of 300 kfci as a function of 
tip height zq. The modeling results in this Figure are labeled by the length L\ of 
tip that is allowed to reorient its magnetic moment. There are three parameters 
in this analysis - a global multiplicative factor, M t M m / K\, and the reduced domain 
wall energy a w j \±§L\K\. Figure Hh plots the best fit values for M t M m /Ki and 
X 2 = T,{^ experimental- AS mode i) 2 / (N-l) (N the number of data points). In all cases the 
best fit value for o w j \i§L\K\ is 0, and the x 2 value at o w j \i§L\K\ = 1 is approximately 
double that when a w j ' \i§L\K\ = 0. Increasing the domain wall energy requires larger 
switching fields: the best fit value at L\— 64 nm for M t M m / K\ increases from 20.8 to 
28.6 when a w j \±§L\K\ increases from to 1. The domain wall energy for Co is reported 
to be a w = 25 ± 3J/m 2 [8]. This leads to a w j \i§L\K\ = 0.85, using L\ = 48nm and 
Xi=0.25 M t 2 (neglecting crystalline anisotropy), with M t = 1.4 x 10 6 A/m, so that our 
fits are consistent with the calculated wall energy, if one allows for a doubling of the best 
X 2 value. The tip end may be magnetically poorly coupled to the rest of the tip because 
of an inhomogeneity or grain boundary: it appears (Figured] inset) to have granularity 
on the scale of a few tens of nm and a kink about 50 nm from its end. We have observed 
qualitatively similar spatial frequency doubling using several nanotube tips, particularly 
in the smallest diameter tips, where inhomogeneities and weak magnetic coupling are 
fundamentally more difficult to avoid. 

The shift Auu in the resonance frequency ujq of the cantilever due to a force gradient 
between the tip and sample dF z /dz is given by A(uj)/uj = —{dF z /dz^)/2k ) where k 
is the spring constant of the cantilever. The phase shift S of the cantilever response is 
then given by 

tan5 = — i—r, (10) 

where Q is the quality factor, uj is the driving frequency and uj' is the perturbed 
resonance frequency of the cantilever. At z = 15nm our model predicts an excursion of 
A(dF/dz)//LLoM t M m w 2 x 10 _9 m. Using a driving frequency at optimal sensitivity 
uj = v (l - 1/V8Q), k = 2.8N/m, Q=350, uj = 168kHz, and estimating M t = 
1.4 x 10 6 ^4/m [9] and M m = 2 x 10 5 A/m [10], this corresponds to an excursion 
in the phase shift AS = 0.7°, in reasonable agreement with the experimental value 
of AS = 0.4° given the uncertainties in the values for the saturation magnetizations. 
Using M t \ < H x > c \/K x <~ 2 [5j, the best fit value M t M m /Xi - 17 (Figure ED) 
implies a critical field of approximately 2.4xl0 4 A/m, much smaller than the saturation 
magnetization of cobalt of 1.4xl0 6 A/m, but comparable to a switching field of 3.2xl0 4 
A/m reported for 30 nm thick, 0.34 /im wide, 2.04 /im long ellipsoidal amorphous 
cobalt nanodots pj]. This reduction in switching field could result from competition 
between the crystalline and shape anisotropics [9] if, for example, the uniaxial crystalline 
anisotropy favors moment alignment along the tip radial direction, while the shape 
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Figure 4. (color online) (a) Full-scale variation in phase angle along cross-sections 
through the centers of the recorded tracks in Figure [2)3. The + symbols represent 
experiment. The other symbols represent modeling as described in the text. The 
modeling curves are labeled by the length L\ of tip that switches magnetic moment 
orientation, (b) Best fit values for M t M m /K\, and \ 2 = Y2(A5 exp — AS mo d) 2 / \N — 1) 
as a function of the tip switching length Li, with a w j ^L\K\ = 0. 



anisotropy favors the axial direction. The tip material could also have a complicated 
structure incorporating grain and domain boundaries, reducing the anisotropy energy. 
The highly non-uniform fields in our case could also play a role in reducing the switching 
field. 

Although our analysis has been presented using a specific model for the tip 
dynamics, the conclusion that the tip moment flips at relatively small fields can be 
presented simply: The magnetic field at the tip must be less than the magnetic field 
at the medium surface, which is given by the saturation magnetization of the medium. 
The fields required to flip the tip are expected, using for example the Stoner-Wohlfarth 
model [5], to be about the saturation magnetization of the tip, which is for epitaxial 
cobalt much larger than the saturation magnetization of the medium. The switching 
fields of our nanotube, just as for amorphous nanodots, are much smaller than those 
of epitaxial Co nanodots [3, [9], which are comparable to the saturation magnetization 
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of cobalt. It might be possible to avoid switching and the attendant spatial frequency 
doubling by developing processes to epitaxially coat the nanotube or by using single- 
crystal nanorods. However, such tips would also generate larger local magnetic fields 
at the sample, increasing the possibility of changing the magnetic state of the sample, 
particularly for the smallest samples. Our results show that reliable information on the 
moment orientations of the media can be inferred even in the presence of tip switching 
if it is realised that the MFM signal is proportional to the absolute magnitude of the 
tip-sample force gradient. 
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